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A complete adiabatic transport of Bose-Einstein condensate in a double-well trap is investigated 
within the Landau-Zener (LZ) and Gaussian Landau-Zener (GLZ) schemes for the case of a small 
nonlinearity, when the atomic interaction is weaker than the coupling. The schemes use the constant 
(LZ) and time-dependent Gaussian (GLZ) couplings. The mean field calculations show that LZ and 
GLZ suggest essentially different transport dynamics. Significant deviations from the case of a 
strong coupling are discussed. 
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I. INTRODUCTION 



Nowadays the trapped Bose-Einstein condensate 
(BEG) is one of the most active topics in modern physics, 
see monograph fT] and reviews [l, H, 0, [E Hi • Between 
many branches of this activity, investigation of weakly 
bound condensates or multicomponent BEG is of keen 
interest. It involves various aspects of Bose Josephson 
junction m, ^ in cluding those in periodically modulated 
traps S IMlni and optical lattices Jl [H [l3], trans- 
port problems in double- well [fO, El) Efi | , triple- well 
[ll[B, [il[ii,[23j23j2l|2|, and multi-well [l| traps, 
topological states [2^ [2a] , etc. In all these areas the non- 
linearity caused by interaction between BEG atoms plays 
an essential role. It drastically enriches BEG dynamics 
by new effects and phenomena. At the same time, the 
nonlinearity complicates and even hampers some process, 
e.g. the adiabatic transport [Tsl. [l9l[2Q| . 

In this paper, we investigate the influence of a week 
nonlinearity on the adiabatic transport in a double-well 
trap. The transport is produced by a controlable and 
irreversible tunneling through the barrier separating the 
wells. It assumes that BEG atoms, being initially in one 
potential well, are completely transferred to another well 
and then kept there. The process can be driven by vary- 
ing the system parameters in time, such as space sep- 
aration between the wells and relative position of the 
well depths. Such transport can be realized in in multi- 
well traps [1] and arrays of selectively addressable traps 
[27| . In this problem we actually deal with two weakly 
bound (through the barrier) condensates or multicompo- 
nent BEG with the components defined as the popula- 
tions of the wells. Being produced, the complete irre- 
versible transport could serve as a useful tool for general 
manipulations of the condensate. Besides, it could open 
interesting perspectives for generation and investigation 
of various geometric phases |20l l28j , creation of topolog- 
ical states Ip5ll26|. etc. 

The transport in multi-well traps can be produced 



by many ways: Landau-Zener [ij, |l5|, ll7| and Rosen- 
Zener [2^ methods, periodic time-dependent potential 
modulation [l^l, Rabi switch |23 | . Stimulated Raman 
Adiabatic Passage (STIRAP) [M [11, H [SI, etc.. 
All these methods provide a robust population transfer 
of the ideal (without interaction) condensate but often 
suffer from the detrimental influence of the nonlinear- 
ity [H, [ii, HO, m HI]. This especially concerns the 
adiabatic population transfer methods, like STIRAP for 
triple- well/level systems [l^, HO, HH , which, being gener- 
ally robust to small variations of the process parameters, 
are, nevertheless, fragile to the nonlinearity. Then a nat- 
ural question arises: is it possible to turn the nonlinearity 
from the detrimental to favorable factor of the adiabatic 
transport? 



As was recently shown [15|, the nonlinearity can in 



deed favor the adiabatic transport if it is modeled within 
the Landau-Zener (LZ) [s^ and Gaussian Landau-Zencr 
(GLZ) [ill protocols. The latter protocol assumes a 
Gaussian time-dependent monitoring of the coupling be- 
tween the wells (barrier penetrability) fl{t) and a linear 
evolution in time of the difference A{t) between the well 
depths. The GLZ is thus a generalization of both LZ 
[33 | and RZ [sj methods. It was shown [ll| that in 
LZ and GLZ the transport is asymmetric, i.e. can be 
considerably enhanced or suppressed by the nonlinear- 
ity, depending on its sign (i.e. repulsive or attractive 
character of the interaction). Similar results were ob- 
tained earlier for LZ-controUed interband transitions in 
accelerated optical lattices [1, [Tsj. This useful feature of 
LZ-based schemes was used in [l^ to build the success- 
ful transport protocols under a strong nonlinearity. It 
was shown that, by choosing the proper monitoring, one 
may produce a complete irreversible transport in wide re- 
gions of repulsive and attractive interaction. Moreover, 
the transport in both left-right and right-left directions is 
possible. Hence the scheme is indeed universal. Though 
LZ and GLZ have much in common, the latter was shown 
to be more flexible. Besides, it suggests a new transport 



regime where not the center but edges of the Gaussian 
coupHng play a decisive role. The LZ and GLZ protocols 
were analyzed in terms of nonlinear structures arising in 
the stationary spectra. 

In this paper we continue the development of the LZ 
and GLZ transport schemes but now for a case of a weak 
nonlinearity UN < fl where N is the total number of 
BEG atoms and U is the atomic interaction. This case is 
not covered in [l^ and hence devotes an additional anal- 
ysis. We will show that, at small nonlinearity, a complete 
population transfer is possible for both LZ and GLZ but, 
unlike the strong nonlinearity case, it takes place at suf- 
ficiently low process rates (which was actually expected 
following the previous studies [H [H [3, 13, El)- At 
even lower rates, i.e. close to the adiabatic limit, the LZ 
transfer is heavily spoiled by strong and slowly damped 
Rabi oscillations while the GLZ protocol is reduced to the 
Rabi switch. The nonlinear structures (loops, etc) can 
appear in the stationary spectra but, unlike the strong 
coupling case [l^, they cannot be already used as re- 
liable tools for the analysis of transport features. The 
transport asymmetry weakens for both LZ and GLZ. 

This paper is outlined as follows. In Sec. II the pop- 
ulation transfer methods are sketched. In Sec. Ill the 
relevant mean field formalism is done. The numerical re- 
sults are discussed in Sec. IV. The conclusions are given 
in the Sec. V. 



II. TRANSPORT PROTOCOLS 

The LZ tunneling between energy levels [s^ is a gen- 
eral physical process which can be straightforwardly re- 
cast for the nonlinear transport of BEG atoms in a 
double- well trap The resulting scheme 

is illustrated in Figure 1 a)-b). 

As seen from Fig. la)-b), the LZ transfer is controlled 
by the constant coupling Q and time-dependent differ- 
ence between the well depths A{t) = at where a is the de- 
tuning rate. Following Fig.lc), in the GLZ protocol, the 
coupling is defined as a time-dependent Gaussian pulse 

nit)=Km, m^eM^^^^}, (1) 

where K is the amplitude, F is the pulse width, and t is 
the centroid time. In the experimental setup, the detun- 
ing can be monitored by varying the well depths while 
the Gaussian coupling by the proper change of the sepa- 
ration distance between the wells [27| . 

In the linear (without nonlinearity) case, the final 
probability of the LZ transfer reads [33] 

P = 1 - e-"^ . (2) 

It allows a complete transition P — I only in the adia- 
batic limit a ^ 0. However, as is seen from the rough 
estimation [13] 

P«l_e-2-ft^(i+4) , (3) 




FIG. 1: a) The general population transfer scheme in a 
double-well trap. The BEG atoms, being initially in the left 
well (bold dot) are then completely transferred in the direc- 
tion indicated by the upper arrow. Ei{t) and E2{t) are the 
ground state energies of the separated wells. The transfer is 
driven by the detuning A{t) and coupling (barrier penetrab- 
ulity) f2. In both LZ and GLZ, the detuning linearly depends 
on time, b) The LZ protocol with the constant coupling Q. 
c) The GLZ protocol with the Gaussian time-dependent cou- 
pling i}{t). 

inclusion of the nonlinearity factor A ~ f7 makes possible 
the complete transport in much wider range of a. More- 
over, the effect is obviously asymmetric with respect to 
the A (interaction) sign. This peculiarity of the nonlinear 
LZ scheme was used in [isj to turn the nonlinearity from 
the detrimental to helpful factor for the transport. 

The introduction of the GLZ protocol is motivated by 
the well-known fact that the LZ transfer actually takes 
place only within a finite time interval near the symmetry 
point A(t) « when A{t) < fi. Then it is natural to use 
a time-dependent coupling of a certain duration, say of 
the Gaussian form. Note that in fact the GLZ protocol 
is a generalization of LZ (constant coupling and time- 
dependent detuning) [s^ and RZ (constant detuning and 
time-dependent couphng) Q schemes. 



III. MEAN FIELD MODEL 

BEG transport is studied in mean-field approximation 
within the Gross-Pitaevskii equation (GPE) [s^ 

ih^ir,t) = [_|iv2 + ^ext(r,t)+5o]«'(r*,i)]']*(r',t) (4) 

where the dot means time derivative, ^'(^^ t) is the order 
parameter of the system, I4xt i^j t) is the external trap 
potential involving both (generally time-dependent) con- 
finement and coupling, go = 4,iTa/m is the parameter of 
the interaction between BEG atoms, a is the scattering 
length, and m is the atomic mass. 

In the two-mode approximation [35| , the order param- 
eter in a double-well trap can be written as 

^'{f,t) = ^{4>i{t)^i{f) + Mt)^2{r)) (5) 
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where <&fc(r) is the static ground state solution of (|4]) for 
the isolated k-th well and 



and 



(6) 



is the amplitude, related to the corresponding population 
Nk{t) and phase (pkit). The total particle number N is 

conserved, i.e. / dr1*(r, t)|V^ = EfcLi Nk{t) = 1. 

Note that the approximation ([5]) is generally valid for 
a weak interaction and small number of atoms, say N < 
1000, see discussion [1, HI]. In this connection, using the 
GPE within the two-mode approximation is somewhat 
contradictory, since the former assumes a large number 
of atoms N while the latter is valid for small condensates. 
Nevertheless, a reasonable balance between these condi- 
tions is possible and their combination is widely used in 
studies of BEG dynamics in double-well traps, see e.g. 

By using the linear canonical transformation [20j 



z^Ni-N2, Z = Ni + N2 



(7) 



(15) 



However, in the present study they come as input pa- 
rameters. In GLZ, the coupling is approximated by the 
Gaussian function ([1]). The energies Ei^2 enter the detun- 
ing pT|) and the interaction U between the atoms inside 
every well is included into the nonlinearity parameter 

m- 

In the previous study [20|, the stationary spectra 
(chemical potentials) 

/i = i[A(t)z + Az^ -n{t)y/l - z2 cos2e'], (16) 

were employed for the analysis of the transport. These 
spectra are obtained by using the equations z = 9 = 
which explicitly read 



(17) 



(8) 



it is convenient to turn the unknowns Nk and 0^ to new 
variables, population imbalance z and phase difference 9. 
This allow to extract from the equations the integral of 
motion Z and corresponding total phase O. 

Then, by substituting dS])-® into performing the 
spatial integration, and using ([T])-®, we obtain equa- 
tions of motion 0, [20| 



z'- sm 



^ = i[A(i) + Az-ff^(t)-^ 
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: cos 26*1 



(9) 



(10) 



where i7(t) — Vl{t)/K is the normalized coupling ((T|), 



m = l^{E,{t) - E2{t)) ^ at 



is the scaled detuning, 



A = 



UN 
'2K 



(11) 



(12) 



is the key nonlinearity parameter determining the ratio 
between the coupling amplitude K and interaction U . In 
©-(Uni), the time t is rescaled as 2Kt — > t. Eqs. @- 
pO)) are solved to get the numerical results presented in 
the next section. 

In principle, the values l^(i), Ei 2{t) and — U are 
determined from the GPE as 

n{t) = -^Jdf [|^v$j • v$2 + ^v^^timi] , (13) 



(14) 



A{t) + ziA+{-iy 



r) = o 



(18) 



where n is an integer real number. Substituting numeric 
solutions of p?l) - p^ into p^ . we get the chemical po- 
tentials fi- and /x+ as the eigenvalues of the stationary 
states. They will be used below to demonstrate the non- 
linear structures arising in the stationary spectra. 



IV. RESULTS AND DISCUSSIONS 

Results of our calculations are presented in the Figs. 2 
and 3. In Fig. 2, the populations of the first well, Ni{t), 
and the second well, N2{t), are shown for LZ and GLZ 
protocols. BEG atoms are initially placed in the first well, 
e. g. Ni{t —00) = 1, N2{t — > —00) — 0, and then 
transferred to the second well. The weak nonlinearity 
A = ±0.5 covering both repulsive (A > 0) a.nd attractive 
(A < 0) interaction is used. Following [isl, HO] , a weak 
nonlinearity allows a complete adiabatic transport only 
for rather small detuning rates a. Here we consider even 
smaller rates 0.01 < a < 0.1 to test the vicinity of the 
adiabatic limit. As discussed below, just this rate region 
demonstrates the major differences between LZ and GLZ 
protocols and is thus most interesting for our aims. Fol- 
lowing Fig. 2, these differences vanish while approaching 
the upper value a — 0.1, where the complete transport 
is most robust, see panels i)-l). At even higher a, the 
process is already too rapid to keep the adiabaticity at 
the given A and the transport becomes incomplete [20| . 

Figure 2 shows that the major LZ-GLZ differences take 
place at the lowest rate a — 0.001 (panels a)-d)). The LZ 
demonstrates here high-amplitude slowly-damped Rabi 
oscillations of the populations. This is because, at so 
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FIG. 2: Populations A''i (dash line) and (solid line) at slow detuning rates a and weak nonlinearity A = ±0.5 for LZ (left) 
and GLZ (right) protocols. Positive and negative A represent repulsive and attraction interactions, respectively. 



small rate, the energies Ei{t) and E2{t) are almost par- 
allel and the condition A < for Rabi oscillations is 
fulfilled for very long time. The transport process in this 
case is vague. However, as seen from panels e) and g), it 
becomes more distinctive with increasing a. 

Much better transport results take place for GLZ (pan- 
els b),d),h)) where the Rabi oscillations exist only for a 
short time determined by the duration of the Gaussian 
coupling fl{t). In this case, the completeness of the trans- 
port depends on the instant when the coupling is over 
and Rabi oscillations are switched off. In other words, it 
is determined by the width T of the Gaussian coupling 
pulse Panels b), d) and d) show that if the Rabi 
switch is done at the proper time, then the transport is 
effective. If not (panel f)), then the transport fails. The 
process can be controlled by monitoring both T and a. 
It is similar to other Rabi switch techniques, e.g. [13]. 

Panels i)-l) show that at higher rates, namely at a = 
0.1, both LZ and GLZ demonstrate a robust and com- 
plete transport (with some advantage of the GLZ proto- 
col). The convergence of LZ and GLZ results is explained 
by the fact that at this rate the condition A < is kept 



in LZ already for much shorter time comparable with 
the GLZ Gaussian pulse duration. This time interval is 
already not enough for development of high-amplitude 
Rabi oscillations. The oscillations also vanish in GLZ 
thus signifying the conversion of the Rabi switch mecha- 
nism to the robust adiabatic population transfer. As was 
mentioned above, a considerable further increase of a is 
not desirable since then the process will be too rapid to 
support the adiabatic following. 

Comparison of the results in Fig. 2 for the repulsive 
(A — -fO.5) and attractive (A = —0.5) interaction shows 
that the asymmetry in BEG transport, i.e. detrimental 
or favorable effect of the nonlinearity depending on the 
interaction sign, is generally small (though it can mani- 
fest itself in particular fragile cases of the Rabbi switch, 
see panels f) and h)). At least the effect is much weaker 
than for a strong nonlinearity when it causes a drastic 
support or suppression of the transport [20| . Since the 
asymmetry effect depends on the nonlinearity magnitude, 
its insignificance for a weak nonlinearity is natural. 

Altogether, Fig. 2 allows to conclude that the cases of 
weak and strong nonlinearity are quite different. Though 



5 



0,8 
0,0 



Constant coupling 







0,6 
0,0 
-0,6 



c) 



-70 



0,5 
0,0 
-0,5 



-40 



0,5 
0,0 
-0,5 



g) 



-40 

0,8i - 



0,0 

-o,a 



10 



0,6 
0,0 
-0,6 



-10 



70 



70 



40 



40 



10 



10 



A=0.5 

a=0.01 



A=-0.5 
a=0.01 

A=0.5 
a=0.025 

A=-0.5 
a=0.025 

A=0.5 
a=0.1 



A=-0.5 
a=0.1 



Gaussian coupling 




-0^65]- 
-40 
0,65i 



0,00 



-0,651- 



-40 

0,8| 

0,0 j) 



-0,8i- 



-10 



0,6i 

o,o' 



-0,6 



-10 



40 



10 



10 



FIG. 3: Chemical potentials /U- (upper lines and structures) and /tt+ (lower lines and structures) for LZ (left) and GLZ (right) 
protocols. The /i_ is depicted by solid (dash) curves for repulsive (attractive) BEG and vise versa for The detuning rate 
a and nonlinearity A have the same values as in Fig. 2. 



the complete adiabatic transport is possible in both cases, 
a weak nonlinearity is distinguished by an essential role 
of the Rabi oscillations (hence specific transport regimes 
like the Rabi switch) and faint asymmetry effect. 

The peculiarities of weak nonlinearity are additionally 

demonstrated in Fig. 3 where stationary spectra /i+ and 
II- are depicted. It is seen that the GLZ spectra exhibit 
nonlinear structures (panels b),d),f),h)) despite a weak 
nonlinearity. However, these stationary spectra in gen- 
eral and nonlinear structures in particular do not display 
the crucial role of Rabi oscillations pertinent to a weak 
nonlinearity and so can hardly be used for the reliable 
treatment of the transport regimes. In this sense, the 
LZ stationary spectra are not instructive as well. They 
do not exhibit nonlinear structures at all and so might 
assume a robust transport. However, as seen from Fig. 2 
a),c),c),d), g), the LZ transport is heavily damaged by 
the Rabi oscillations. Here we see again a big difference 
with the case of a strong nonlinearity where the solid 
nonlinear structures in the stationary spectra allow to 
do a reliable analysis of the transport features. 



V. CONCLUSIONS 

The complete adiabatic transport of Bose-Einstein 
condensate in a double-well trap is investigated within 
the conventional Landau-Zener (LZ) protocol and its gen- 
eralization to timc-depcndcuit Gaussian coupling (GLZ) 
for a case of a small nonlinearity. The relevant range 
of slow detuning rates when such nonlinearity manifests 
its main peculiarities is analyzed in detail. The essential 
role of Rabi oscillations in the vicinity of the adiabatic 
limit is demonstrated. These oscillations hamper the LZ 
transport but make possible the Rabi switch transport 
for GLZ. For higher detuning rates, both LZ and GLZ 
converge to a robust adiabatic transfer. The asymmetry 
effect is found generally small. 

The nonlinear structures in the stationary spectra arise 
at very low detuning rates despite a small interaction. 
However, these structures do not reveal the significant 
role of Rabi oscillations and so, unlike the case of a strong 
nonlinearity, cannot be used as a reliable tool for the 
analysis of transport features. 
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Altogether, the calculations show a considerable differ- 
ence between slightly and strongly nonlinear transports. 
They demonstrate different mechanisms. What is most 
important, a weak nonlinearity with its negligible asym- 
metry effect cannot be used as a powerful tool to enforce 
or suppress the adiabatic transport. In this sense a strong 



nonlinearity is certainly more promising [20|. 
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